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Abstract 

We study regular, static, spherically symmetric solutions of Yang-Mills 
theories employing higher order invariants of the field strength coupled 
to gravity in d dimensions. We consider models with only two such in- 
variants characterised by integers p and q. These models depend on one 
dimensionless parameter a leading to one-parameter families of regular so- 
lutions, obtainable by numerical solution of the corresponding boundary 
value problem. Much emphasis is put on an analytical understanding of 
the numerical results. 
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1 Introduction 



Regular, static, and finite energy solutions to the Einstein- Yang-Mills (EYM) 
equations have presented a subject of interest since the work of Bartnik and 
McKinnon [T] (BM) in 4 spacetime dimensions. Besides the regular solutions 
also corresponding black hole solutions were found [3 El El representing novel 
examples of hairy black holes. In addition, these classical solutions are of interest 
also as the 'solitons' of the low energy effective actions of superstring theory [Hj, 
since both gravitational and non Abelian gauge fields occur in the latter and 
in gauged supergravities. More specifically, such solutions (especially the black 
holes) have an important role to play in quantum gravity and also on (higher 
dimensional) .D-branes [7j. In this wider context therefore, the construction of 
such solutions in spacetimes of dimensionality higher than d = 4 is of actual 
interest. This is the purpose of the present work, namely the study of EYM 
systems in arbitrary spacetime dimensions d. 

Although the YM theory in flat space has no 'soliton' type solutions due to 
its scaling properties, the coupling to gravity changes the picture. Bartnik and 
McKinnon discovered a discrete family of static, spherically symmetric globally 
regular solutions. However, the scaling properties of the EYM theory expressed 
by a 'Virial Theorem' forbid such solutions for d > 5; while the scaling dimension 
of static gravitational fields is d — 3 the one of the YM fields is equal to d — 5. 
The virial theorem implies that non-trivial regular solutions require terms in the 
action with scaling dimensions of different sign. The way out proposed in El is 
to use higher order invariants for gravity and/or YM fields. The choice of models 
boils down, in principle, to the selection of various terms in the gravitational 
and non Abelian hierarchies of increasing orders in the respective curvatures, 
namely the Riemann and the YM curvatures, which are reparametrisation and 
gauge invariant. There is an important restriction, namely that one considers only 
those Lagrange densities that are constructed from antisymmetrised 2p curvature 
forms, and exclude all other powers of both Riemann and YM curvature 2-forms. 
(In the gravitational case this results in the familiar Gauss-Bonnet Lagrangean.) 
As a result, only velocity-squared fields appear in the Lagrangean, which is what 
is needed for physical reasons. Here we will add only the minimal number of 
such higher order terms that are necessitated by the requirements of the virial 
theorem. This criterion makes the inclusion of higher order gravitational terms 
unnecessary and in fact the inclusion of higher order gravitational terms does 
not seem to alter the qualitative properties of the classical solutions From 
the viewpoint of classical solutions, indeed, such terms seem redundant. They 
become, however, important in the mathematically interesting situations where 
the first member of the YM hierarchy, namely the usual YM term, is absent. 
The situation concerning higher order YM curvature terms in the string theory 
effective action is complex and as yet not fully resolved. While YM terms up 
to F A arise from (the non Abelian version of) the Born-Infeld action [TJIj, it 
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appears that this approach does not yield all the F 6 terms [TT]. Terms of order 
F 6 and higher can also be obtained by employing the constraints of (maximal) 
supersymmetry ^2]. The results of the various approaches are not identical. 

Thus we end up adding only higher order YM terms to the usual EYM La- 
grangean. The YM hierarchy employed is 

F(2p) denoting the p fold totally antisymmetrised products 

F(2p) = = F A F A . . . A F , p times , (2) 

of the YM curvature, F(2) = F^ u , in this notation. Clearly, the highest value P 
oip in Eq. (JJJ is finite and depends on the dimensionality d of the spacetime. To 
complete the definition of the models (JTJ) the gauge group G must be specified. 
With our aim in the present paper, of constructing static spherically symmetric 
solutions in d spacetime dimensions, the smallest possible such gauge group is 
G = SO(d — 1) requiring d > 4. 

In (HI IH] the simplest possibility P = 2 was chosen. The scaling dimension of 
the p = 2 term is equal to d — 9, yielding a negative contribution in the virial 
theorem for d < 9. As a consequence these gravitating YM models in spacetime 
dimensions d — 5, 6, 7, and 8 were found to support particle like solutions 1 . Here 
we will consider more generally models involving two arbitrary terms of the sum in 
Eq. (£[}, say F(2p) and F(2q). For simplicity let us call them p-g-models. Like the 
1-2-model these simplest nontrivial models possess one dimensionless parameter 
a providing a measure of the strength of the gravitational self-interaction of the 
YM configurations. 

The regular solutions reported in [HI E] , which were constructed numerically, 
featured numerous novel and exotic properties, some of which were not under- 
stood or explained at a more fundamental level. The purpose of the present work 
is specifically to understand these properties analytically, and then substantiate 
the conclusions and predictions numerically. In fact, our analysis immediately 
generalises to the general p-g-models, all showing the same basic structure. In 
this context it also turns out to be useful to interpolate in the dimension, i.e. take 
the dimension parameter d to be an arbitrary real number. This is possible since 
the reduced one dimensional equations of motion for static, spherically symmetric 
configurations depend parametrically on d. 

The properties of the solutions alluded to in E] are: 

1. The solutions form one parameter families parametrised by a. It was found 
that a increases from the flat limit zero up to some a max and then turns 
back. 

1 Taking only the p = 2 term one finds generalisations of the BM solutions in these dimensions 
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2. In d = 6, 7, and 8, where there is a finite-energy gravity-decoupling solution, 
two solutions appear to exist for most values of a, up to a max . Starting 
from the flat F(2) + F(4) solution a increases from zero to a maximal value 
ct max and then runs back to zero for d — 8, while for d = 6 and 7 the families 
end at some finite non-zero value of a 

3. In d = 5 there is no flat F(2) + F{4) solution, but the relevant flat solution 
is the instanton of the pure F(2) theory considered as a static solution of 
the d — 5 theory After a suitable rescaling both gravity and the F(4) 
term decouple in the limit a — > 0. Again a turns back after reaching some 
maximal value, but then seems to converge to some limiting critical value 
a c after performing a finite or infinite number of oscillations. 

4. Unlike in the case of d — 4 EYM theory no multi-node solutions were found. 

In the present paper we will try to clarify the findings desribed above with 
particular emphasis on the critical values of the parameter a marking the end- 
point of the one-parameter families of solutions. Based on experience with other 
systems like self-gravitating magnetic monopoles ^Hj we propose that the critical 
solutions run into certain fixed points (f.p.) of the dynamical system represented 
by the reduced field equations, a system of nonlinear ordinary differential equa- 
tions. Employing the standard method of linearisation at the f.p. it is possible 
to study the analytical properties of solutions coming close to the f.p. and derive 
asymptotic expressions for their parameters. Constructing solutions starting at 
the f.p. itself provides a direct method to determine critical parameters like a c 
yielding much more precise values than extrapolation from solutions coming only 
close to it. In fact, our analysis leads to a rather clear picture of the results listed 
above. 

Concerning point 2., (d = 6, 7, and 8) we claim that the upper branch of the 
solution family runs back all the way to a = 0, the limit (after suitable rescaling 
decoupling the F(2) part) being the analogue of the BM solution for the F{4) 
model. As a tends back to zero the solutions for d = 6 and d = 7 come very 
close to a fixed point already known from the EYM theory in d = 4 related to the 
behaviour at the horizon of the extremal Reissner-Nordstr0m black hole [14^ I13j. 
Our analytical analysis, however, excludes a bifurcation with solutions running 
into the fixed point. Although our improved numerical computations support 
this analytical result for d = 7, for d = 6 the approach to the f.p. is too close to 
resolve this question numerically. 

The limiting solution for d — 5 described in 3. is related to a new type of fixed 
point described in Section El The resulting structure is comparable to the limit 
of infinite Higgs mass for magnetic monopoles [T3]. 

As to 4. we show that multi-node solutions exist, but only for d between 4 
and 5, i.e. not for the integer dimensions studied in [HI El • 
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In addition to clarifying the points 1. — 4. above, we have carried our conclu- 
sions further, to verify that these qualitative features repeat themselves for the 
general p-g-models. Particularly the models with q = p + 1 show a certain degree 
of periodicity modulo 4p + 1. 



2 The reduced one dimensional theory 

The Lagrangean of the Yang-Mills hierarchy in d spacetime dimensions intro- 
duced in jS] including the generic p-th member, is 



YM 



E 



^2.(2p)! 



trF(2p) 5 



(3) 



where the 2p form F(2p) is the totally antisymmetrised product of the YM field 
strength F(2). For the metric the static spherically symmetric ansatz with d — 1 
space-like dimensions is 

ds 2 = ds\ - r 2 dn 2 d _ 2 , (4) 

where ds\ is the metric on the two dimensional orbit space factorising out the 
action of the rotation group and dVl? d _ 2 the invariant line element of S d ~ 2 . The 
2d metric can always be brought to the diagonal form 

e 2 »dt 2 - e 2X dR 2 . 



dsl 



(5) 



Considering only time independent solutions we naturally take t to be the Killing 
time, while the radial coordinate R remains arbitrary. Plugging this ansatz into 
the standard Einstein-Hilbert action results in the reduced one dimensional La- 
grangean 



v+X „d-A 



d-3 + e- 2A (i/(r 2 )' + {d-3)(r'f 



(6) 



where G denotes Newton's constant. The static spherically symmetric Ansatz for 
the SO±(d) resp. SO±(d — 1) in even resp. odd spacetime dimensional YM field 
is 

i/i±iV 



A n = 



A; 



1 - w 



where T c h is the chiral matrix in the appropriate dimensions. 
The resulting reduced one-dimensional Lagrangean is 



E r 

V 



d-4 



yp-i 



a p e v - x (W') 2 + b p e u+x r 2 V 



[r„r,], (7) 



(8) 



with V = (W 2 — l) 2 /r 4 and coefficients 
(rf-2)! _1 A 

° p ~ Tp 2(2p-iy.(d-2 P -iy. ~ 2 Tp ' 



and 



d-2p 



2p 



(9) 



From the form of these cofficients one can see that if d is an integer it must be 
larger than 2p for the term with F(2p) to contribute. Thus e.g. for d = A only the 
F(2) term is present. Therefore we will restrict our analysis to values of d > 2P, 
where P is the maximal value of p in the sum. 

An infinitesimal scale transformation of the action S = J dR^La+LyM) yields 
a virial theorem. While the gravitational part picks up a factor d — 3 the F(2p) 
terms get a factor d — Ap — 1. As a consequence non-trivial flat solutions exist 
only, if the conditions d — AP — 1 < and simultaneously d — Ap — 1 > for 
the first non-vanishing term in the YM part are fulfilled. On the other hand, 
if d — Ap — 1 < for the first non- vanishing term then the gravitational part is 
required to obtain non-trivial solutions. 

Varying the action with respect to the variables v, r and W yields three 
second order ordinary differential equations. Introducing independent variables 
N, k and U for the first derivatives 

N = e _ V , k = reals' + N , U = e- x W , (10) 

the equations of motion become 

re- x N' = (K-N)N -2Gj2a P V p ~ 1 U 2 , (11a) 

p 

re~V = — (d - 3)(d - 6) - k 2 + -(d - 4)(d - 5)iV 2 (lib) 

_ G Y^V p - l {a p {d-Ap-2)U 2 -b p {d + Ap-$)T 2 ) , (11c) 



re~ x U' 



p 
1 



Y J V p ' 1 {{Ap+l-d)N -n)a p U (lid) 

(lie) 



E p a p VP- 1 

+2WT £ V p ~ l (pb p - (p - l)a p ^j 



2 

I 



where we have introduced the shorthand T = (W 2 — l)/r. The variation with 
respect to A yields the algebraic constraint for the first derivatives 

kN --(d-3) + -(d-5)iV 2 -GY.V^UpU 2 -b p T 2 ) = 0. (12) 
2 2 p 

Subsequently we shall restrict ourselves to models with only two terms (pl,p2) 
in the sum over p denoting them simply as p and q. The simplest choice is p = 1 
and q = 2, the choice made in [SI E] • Suitably rescaling r and the total action the 
model depends only on the dimensionless parameter a 4< - 9_p - ) = G 2 ( q - pS) f p 2q ~ 1 f q 1 ~ 2p . 
Taking more terms in the sum does not seem to change the overall picture, but 
makes the numerical as well as the analytical analysis much more cumbersome. 
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3 The F(2p) model 



Before we turn to the F(2p) + F(2q) model we consider the model with only one 
F(2p) term, generalising the system studied by Bartnik and McKinnon [Ij. In 
this case the parameters G and f p can be completely absorbed by a rescaling of 
r and the total action. We still have the freedom of choosing a suitable radial 
coordinate. As in we impose the condition e A = r, call the corresponding 
radial coordinate r and denote r derivatives by a dot. Thus we obtain from 
Eqs. dm) 

r = rN , (13a) 
W = rU } (13b) 
N=(k- N)N - V P ^U 2 , (13c) 

.=i-^+ (rf - 4 f - V-i) 

+ ^((4 P + 2-^ + (d "^" 1 f + 4P " 8 M , d3d) 



2p 

?7=((4p+l-d)JV-K)t/ + iyr(d-2p-l-2(p-l)^) , (13e) 



and the constraints 

- 1 - rT = , (14a) 
2«N + 3 - d + (d - 5)iV 2 - V th - 1 (U 2 - d ~ 2p ~ l T 2 ) = . (14b) 



It turns out to be convenient to introduce new variables y = V p 1 T 2 and 
z = U/T and rewrite the equations as 

f = rN , (15a) 

W = rTz, (15b) 

f=(2Wz - N)T , (15c) 

N=(k-N)N -yz 2 , (16a) 

* = W+ (d - 4 f "V-i) 

+ ( ( 4, + 2-^ + ^-^-^ + ^- 8 ) )|, (16b) 

y=(±pWz- (4p-2)N)y , (16c) 

z=(d-2p-l)W+ ((Ap + 2 - d)N - k)z - 2pW z 2 , (16d) 
and the constraints 

2kN + 3-d+(d-5)N 2 - (z 2 - d ~ 2p ~ 1 )y = 0, (17a) 

v 2p ' 
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W 2 -l-rT = 0, (17b) 
T 2 P _ r 2 P -2 y _ g ^ 17c ^ 

4 Fixed Points of the F(2p) model 

4.1 r = and r = oo 

Looking for globally regular solutions existing for < r < oo we have to cope 
with the singularities of Eqs. (j!3)l at r = and r = oo. Since for regular solutions 
the space time is locally flat at r = and r = oo, we have to require N — k — 1 
there. This implies that the radial variable r behaves locally like lnr and thus 
t — > =poo for r — > resp. r — > oo. With a suitable choice of new dependent 
variables both singular points become hyperbolic fixed points (f.p.s). From this 
viewpoint regular solutions correspond to the 'stable manifold' of the respective 
fixed point [E]. As in the case of the BM solutions regularity requires W — ±1 
and U = at r = and oo. Let us first turn to the point r = 0. Although a 
Taylor series expansion using the 'Schwarzschild' coordinate r as an independent 
variable yields essentially one arbitrary parameter b from W = 1 — br 2 + 0(r 4 ) as 
in the case of the F(2) theory PIE], the linearisation at the f.p. is slightly more 
difficult due to the appearance of the singular expression V = (W 2 — l) 2 /r 4 and 
the auxiliary variable z in the equations. 

For W = N = k = 1 the equation for z has the f.p.s z = 1 and (2p+ 1 — d)/2p. 
It turns out that z = 1 for solutions regular at r = 0, while z = (2p + 1 — d)/2p 
yields a singular mode. This suggests the introduction of the variables 

Ttr W-\ Z-l - N-l , K-l 

W = — , z = , N = , and R = . (18) 

While z, N, and R vanish at r = the variable W has the finite limit —6, where 
b is the parameter of the Taylor expansion. From Eqs. (fT^|) we obtain 

W = 0(r), (19a) 

Z=-dz + (4p + l-d)N -R + 0(r) , (19b) 

ti=-N + R + 0(r) , (19c) 

k={d-4){d-5)N -2R + 0{r) . (19d) 

Replacing the r derivatives by rd/dr = (1 + 0(r))d/dr these equations take 
the form required by Prop. 1 of It is straightforward to see that besides 

the regular r 2 mode for W the linearisation yields a singular one behaving as 
r (2p+i-d)/p_ fj^g (jy^ ft) system yields an eigenvalue 3 — d, and a second one d — 6 
incompatible with the constraint Eq. ()14b ). Thus for d > 3 the behaviour of N 
and k at r = is determined through the non-linear terms in Eqs. (|19|). One 
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finds 

iV=l-^lAr 2 + 0(r 4 ) , (20a) 

K = l + (4p- 3) ^~r 2 + 0(r 4 ) . (20b) 

The linearisation for r — > oo proceeds along similar lines. One finds that 
regular solutions behave like 

W = ±{l- cr (2p+1 - d)/p + ...), (21) 

and there is a singular r 2 mode. Thus the regular and singular modes of W 
exchange their role going from r = to r = oo as to be expected from r — > ±oo. 
Similarly the eigenvalue 3 — d from the (N, k) system now yields a regular mode 

Tfl 

JV = 1 -2^ + -. ( 22 ) 

and solutions regular at r = oo are described by two free parameters, c and the 
'mass' m. 



4.2 W = U = 

Besides the f.p.s for r = and r = oo there is also a f.p. for finite r = tq 
with W = U = N = and « = «o- This f.p. is characteristic for the extremal 
Reissner-Nordstr0m solution at its horizon and therefore we call it the RN f.p. 2 
From the constraints Eqs. (|14|) we find 

^- (d-fr-Vs) ' md = <*-W- 3)- P3) 

Putting r = r (l + 5r), r{7 = f7 and k = k,q + 5k and keeping only linear terms 
we obtain from Eqs. (jllj) the system (eliminating 5r via Eq. ([14))) 

iV = k A^ , (24a) 
6k=-2k q 6k- {d + Ap-8)K N , (24b) 
W = U , (24c) 

j7 = -(d-2p-l)W-Kol7. (24d) 



2 The RN field configuration arising in the d = 4 EYM theory is the infinite node limit 
of the regular solutions, in the external region the function W (r) vanishing everywhere. In 
that case the 50±(4) = SU±(2) gauge group breaks down to SO (2) = 17(1), namely this field 
configuration is an Abelian embedding. In the present cases, the W(r) = configuration of a 
SO±(N) gauge field with N > 5, describes a SO(N — 1) gauge field, namely a non- Abelian 
embedding. 
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The eigenvalues of these linear equations are kq, —2kq for the (N, re) system and 



A± = \ (-K ± ^f(2p-5)d + 2p + 7) , (25) 

for the (W, U) system. The eigenvalues A± are always real for p > 2, while for 
p — 1 and p = 2 they are real only if d < 3 resp. d < 11. 

There is a similar f.p. with iV = re = 1 and W = U = already present for 
the flat F(2p) model. Since then r « lnr this f.p. requires r — * or r — > oo. The 
corresponding eigenvalues A± are 



A- 



- (4p + 1 - d ± ^(d - 4p - 3) 2 - 8p - 4) . (26) 



These eigenvalues are real if d lies outside the interval [4p + 3 — 2y/2p + 1, 4p + 
3 + 2y/2p+ 1]. For p = 2 this interval is [6.528,15.472] and for p = 3 it is 
[9.709,20.292]. As long as the eigenvalues have non- vanishing imaginary part 
the solutions coming close to the f.p. are oscillating. As a consequence we might 
expect to find regular solutions with any number of zeros. However, for d > 4p+l 
the real part of A± becomes negative and thus the solutions coming close to the 
f.p. tend to zero for r — > oo and W cannot reach ±1. In order to analyse these 
oscillating solutions beyond the linear approximation we consider the 'Lyapunov 
Function' 

Z=(z 2 - d - 2 P~ 1 yp^y = (W 2 - lfv- 2 W 2 - d ~ 2p ~\ w 2 - l) 2p , (27) 
v 2p ' 2p 

which is a measure of the amplitude of oscillations (as long as \W\ < 1), and 
satisfies 

Z = 2((4p + 2 - d)N - re) (W 2 - l) 2p ~ 2 W 2 . (28) 

Thus for asymptotically flat solutions with Af— > 1, k — > 1 as r - > oo, the 
amplitude of oscillations decreases for d > 4p + 1, increases for d < 4p + 1, and 
has a finite limit for d — 4p + 1. That means that regular solutions with many 
zeros cease to exist for d > 4p + 1. However, since 4p + 3 — 2y / 2p~+T < 4p + 1 
there is always a finite d interval, where regular solutions with any number of 
zeros are expected to exist. 



4.3 Conical singularity 

Forp > 1 there is another f.p., which turns out to be relevant for the limit b —>■ oo 
in certain dimensions (e.g. d — 5 for the F{2) + F(4) model). It is characterised 
by W = 1 and nontrivial 1) finite values of N, k and z. Putting the r.h.s. of 
Eqs. (fTBjl to zero yields the f.p. equations 

= kN — N 2 — yz 2 , (29a) 
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o=i-^ + (d - 4 f -V-il + ^-y 

+ ^-2p-m +ip s )y 2 (2gb) 

= (Apz - (4p - 2)N)y , (29c) 
= (4p + 2 - d)iVz - KZ + d-2p-l - 2pz 2 , (29d) 

leads to the solution 

Ar2 _ 16p 3 - 8dp 2 + 2pd 2 -8pd + Up + d-3±VQ 

7V °" 2(4p 2 -4p+l)(2p-l) ' W 

k = (2p + 3 - d)iV + 2p(d "^ , (30b) 

(2p- 1)N 

2r> — 1 

^o = -^iVo, (30d) 
ip 

where 

Q = (2pd - 6p + l)(d - 3)(2pd 2 - 16p 2 d - 4pd + d + 32p 3 + lpp - 3) . (31) 

The zeros of the polynomial Q are at d — (Qp — l)/2p, 3, and (16p 2 + Ap — 

1 ± \/12&p'- i — 96p 2 + 16p + l)/4p. In order to obtain a real value for iVg in the 

admissible region the dimension d has to be restricted to the interval 3 < d < d p 

with 

_ 16p 2 + Ap - 1 - V128p 3 - 9Qp 2 + 16p + 1 
^ - — , (32) 

or it has to be larger than (16p 2 + 4p — 1 + y/128p' s — 96p' 2 + lQp + l)/4p. The 
numerical values of d p for p = 2, 3, and 4 are ~ 5.63, 8.63, and 11.81. 

The fact that N has a finite non-vanishing limit iVo < 1 implies that r — > 
for r — ► — oo. As a consequence of iVo 7^ 1 the gravitational field has a 'conical' 
singularity at r = with a solid angle deficit and the YM field has a weak 
singularity 

W = l- \yl /2p r 2 - 1/p + 0(r 2+1 / p ) , (33) 

as for the massive EYM model ('Proca theory') studied in |13j . 

In order to linearise the Eqs. (|TH|) at this 'conical' f.p. we put N = Nq + SN, 
k = no + 6k, y = yo + 8y, z = z + 5z and keep only linear terms in 5N etc. 
Eliminating 5k via the constraint Eq. ([17b) we obtain 

8N=(3- d)N 8N - (| + d ~f ~ 1 )Sy - y z Q 5z , (34a) 
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5'y = -2(2p - l)y SN + 4py 5z , (34b) 
fc = ^((4p - 3)iV + „)M + 3£i(!^fp± - 

+((4-d)JV -K -^ = ^j/o2b)^. (34c) 

Due to the complicated structure of this linear system its eigenvalues can only be 
determined numerically. For small p one finds that there is one pair of complex 
conjugate eigenvalues if d is near d = 2p, but there are three real ones if d is larger 
than some value d^\ For p = 2 resp. p = 3 one finds d% ~ 5.55 and djj ~ 8.0. 
For p 5.2 the value becomes equal to 2p and all three eigenvalues are real 
for all d > 2p. For d = d p one of the three negative modes becomes a zero 
mode, because the two different solutions for Nfi degenerate. As long as there 
are complex eigenvalues the solutions coming close to the f.p. oscillate. 

Since the real parts of all three eigenvalues are negative for d < d p the f.p. 
is always repulsive for r — > — oo. Thus there are no free parameters available 
that would be necessary to construct global solutions regular at r = oo with the 
behaviour ()33j) near r = 0. 



5 F(2p) + F(2q) model; q > p 



Next we study the model with with the two terms F(2p) + F(2q). As mentioned 
above we require d > 2q (and hence d > 2p + 2) and d < Aq + 1. After a suitable 
rescaling this model depends on a single dimensionless parameter a. 
From Eqs. ([lip we get 



r = rN , 
W = rU , 

N=(k-N)N -a 2 V p ~\l + V pq )U 2 , 



(35a) 
(35b) 
(36a) 



1 — K -\- 

2 



(d-4)(d-5) 



(N 2 - 1) 



(4p + 2-d+(4q + 2- d)V pq )U 2 + ( 
(4g-8 + d)(d-2g-l) 



(4p-8 + d)(rf-2p- r 



2p 



2q 



(36b) 



1 



1 + K 



A"/ 



((4p + 1 - d)N - k)U + ((4g + 1 - tZ)N - «)^E7 



ry-2 

+^T(c/ - 2p - 1 + (d - 2q - l)V pq - 2(p - 1 + (q - l)V pq ) — 



J^2 



(36c) 



12 



and the constraint 

2kN + 3 - d + (d - 5)N 2 = a 2 V p ~ l 



(l + V pq )U 2 - ( 
(d-2q-l 



d-2p-\ 
2~p 



2q 



(37) 



with V pq = V q ~ p . 



Again it turns out to be useful to introduce the variables y = V p 1 T 2 and 
z = U/T and rewrite the equations as 



r = rN , 
W = rTz , 

N= ( K - N)N - a 2 (l + V pq )yz 2 , 

2 (d - 4)(d - 5) /Ar2 a 2 

« = i - « 2 + f -(^ - 1) + y 



(38a) 
(38b) 
(39a) 



(4p + 2-d+(4q + 2-d)V pq 



(Ap-8 + d)(d-2p-l) (4q-8 + d){d-2q-l 



+( 



4p(\Vz 



2p 
2p-l 
2p 



+ 



2g 



■V, 



pq 



N y 



(39b) 
(39c) 



4p + 2-d)N - k)z+ ((4:q + 2-d)N -K)V pq i 



1+V p g 

+W(d-2p-l + (d-2q - l)V pq - 2(p + qV % 

with the constraints 

2kN + 3-d+(d- 5)N 2 = a 2 



P'V 



(39d) 



pi j 



(l + ^p 

d - 2p - 1 (d - 2g - 1) 



2p 



+ 



2g 



W 2 - 1 - rT = , 
T 2 P _ r 2p-2 y = q _ 



2/ , (40a) 

(40b) 
(40c) 



The behaviour of solutions regular at r = is now given by the expansions 



N = l-a 2 { 



(4b 2 ) p (4b 2 ) 
4p 4q 

(4b 2 ) p .(4b 2 ^ q - 



)r 2 + 0(r 4 ) 

k = 1 + a 2 ((4p - 3)^- + (4q - 3)^)r 2 + 0(r 4 ) . 



(41a) 
(41b) 



The F(2p) and F(2q) terms in Eqs. (J41j) scale differently. In the limit b — ► oo 
the F(2p) term decouples and one obtains the BM type solutions of the F(2q) 
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2 

model after the rescaling b — > foa 21 ?- 1 and letting a — > 0. Likewise for b — ► 
and d < 4p + 1 one obtains the BM type solutions of the F(2p) model rescaling 

2 

b — ► fret 2 ?- 1 and letting now a — >■ oo. 

From the virial theorem we know that for 4q + l>d>Ap + l there is a 
flat solution, obtainable in the limit a — > 0; the corresponding parameter 6fl at 
diverges as d — > Aq + 1 and vanishes in the limit d — > 4p + 1. 

For ci = 4p + 1 and b — > the F(2q) and gravity terms both decouple and 
one obtains a rescaled flat F(2p) instanton. Linearising the equations for a per- 
turbation of the instanton background yields a finite limit for b 4p ~ 2q ~ l a 2 ; thus a 
vanishes as b — > if 4p < 2q + 1 or diverges if 4p > 2g + 1. 

It is easy to see that for r — > oo the terms coming from F(2q) are subleading 
and hence the asymptotic behaviour is that of the F(2p) model. 

Also the other f.p.s besides r = and r = oo discussed for the F(2p) model 
have counterparts here. Of particular relevance for the application to limiting 
solutions is the RN f.p. with W = U = N = and finite values of r and k 
determined by the equations 

a 2 ( d-2p-\ d-2q-U 
= 3-. + ^(^^ + ^-^), (42a) 



Kl ~ 1 ~ T " /+ 4r 4p " 2 v i> 



(d-4)(d-5) a 2 AAp - 8 + d)(d - 2p - 1) 



+ (49 " 8 ;g-t 29 " 1 ) • < 42b > 



Except for d = 2p + 1 or d = 2q + 1 the resulting equation for r can only be 
solved numerically. 

Similarly there is the f.p. with W = U — 0, N — k — 1 and r — > cxd already 
present in the flat p-g-model. However, for r oo only the F(2p) term survives 
so we can use the corresponding results of the previous section. 

What remains is the conical f.p. of Eqs. (J29|) . It turns out to become relevant 
for the limiting solutions as b — > oo in the interval 2q < d < d q (e.g. d — 5 for the 
1-2- model). The numerical analysis shows that for large values of b the solutions 
regular at r = come very close to the f.p. (of the F(2q) model) within a very 
short r interval of length r$ = l/b q . They consist of an 'interior' part between 
r = and r fa ^/Fq and an 'exterior' one extending from r fa ^/Tq to r = oo (see 
Figs. [Hand El). 

The interior part is dominated by the F(2q) terms in Eqs. (|39J) and the expan- 
sions Eqs. (|41|) suggest the introduction of the rescaled variable f = r/ro = rb q . 
In the limit b —>■ oo the globally regular solutions split into an interior part ex- 
tending from the regular origin f = to the conical f.p. at f = oo and an exterior 
one extending from the f.p. at r = to r = oo. The spacetime corresponding to 
the interior part of the solution has a conical singularity at f = oo with a solid 
angle deficit much like for global monopoles [To] . 
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Solutions for sufficiently large b can be obtained by matching a rescaled inte- 
rior solution to an exterior one in a domain where the linearisation at the f.p. is 
appropriate. Although the F(2p) terms are a small perturbation for the interior 
solution, this is true for the exterior one only near the f.p., while the F(2p) term 
becomes dominant for large r. 

In order to describe the interior part we use the variable f = rb q and introduce 
f = 69 _1 T, such that y = T 2q r 2 ~ 2q = f 2q f 2 ~ 2q . For large b we get from y = 0(1) 
the estimate 

rT 2q-l , „ 

W - 1 « — ~ f— b l ~ 2q . (43) 

Putting Y = (N, y, z) we can linearise Eqs. at the f.p. with Y = Y + 5Yj to 
obtain 

f^SYj = L5Yt + 0(<5y/ 2 ) + 0(f 2 -^b l - 2q ) . (44) 
dr 

L here is the matrix on the r.h.s. of Eqs. (|34jl . with some rescaling due to the 
appearance of a 2 in Eqs. 1)39)1 . In the region 1 <C f 1 / 9 <C 6 we can neglect the 
0(.) terms to get 

(jy 7 = Kf\ = iT(6r) A , (45) 

with some constant matrix K and constant vector rj (A denoting the diagonal 
matrix of eigenvalues). In order that the exterior solution hits the f.p. for r — > 
— oo we have to fine-tune three parameters. Since the solutions regular at r = oo 
depend only on the two parameters c and m introduced in Eqs. (|2T|) and f|22(l. we 
have to tune also a. This determines the 'critical' value of a which in fact is a 
way to compute it more precisely than as a limit of regular solutions starting at 
r = 0. We put £ = (a, m, c) with £ the triple parametrising the solution running 
into the f.p. For £ = £ + 5£ the solution misses the f.p. but comes close to it 
for \8£\ 1. Linearisation at the f.p., valid for r C 1, gives now for the exterior 
solution 

SY E w Kr x M5£ , (46) 

with some constant matrix M. Matching the interior with the exterior solution 
at some fixed r 1 using these approximations gives 

M5i = b qX V . (47) 

If we rescale the S^s with the smallest real part of the A's, the rescaled quantities 
either tend to a limit for b — > oo if the eigenvalues are real, or they oscillate about 
the limiting value in case they are complex. This behaviour is well exhibited in 
Figs. El and H 

6 Numerical results 

The numerical analysis of the F{2) + F(4) model [HI El showed that for d — 5, 6, 7, 
and 8 there are 1-parameter families of globally regular solutions labelled by the 
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parameter a. All these solutions resemble the lowest BM solution with one zero 
of the potential W . In contrast to the BM case no solutions with more than one 
zero were found. All the families have a maximal value a max and end at some 
'critical' value a c , where the solutions cease to exist. While for d = 5, 6, 7 a 
critical value a c > was obtained, it was suggested that a c = for d = 8. Our 
aim is to explain what happens in the limiting cases, determine the values of a c 
and to analyse the corresponding critical solutions. As already mentioned above 
it proves very helpful for this purpose to consider also non-integer values of d, 
interpolating between the physically relevant integer ones. This is made possible 
by the explicit ^-dependence of the Eqs. (ITT]) . In fact we extended our study of 
the F(2) + F{4) model to the whole interval 4 < d < 9 for which globally regular 
solutions can be found. We find that for 4 < d < 5 analogues of the discrete 
family of BM solutions with any number of zeros of W exist, while for 5 < d < 9 
only solutions with one zero were found. 

Let us first discuss the results for 5 < d < 9. Varying the parameter a 
for a fixed d one obtains a smooth 1-parameter family of solutions conveniently 
parametrised by the curves b(d, a) , where b is the parameter characterising the 
solutions at r = 0. These curves start with finite b at a = and depending on d 
either turn back to a = after running through some maximum or end at some 
finite value a c (see Fig. EJ)- Our aim is to characterise the limiting solutions on 
both endpoints of the curve. 

The endpoint b(d, 0) is obtained by switching off the gravitational field taking 
a — > 0. A was discussed in Section El the corresponding solutions are those of the 
flat F(2) + F{4) theory already found some time ago [TZj. Actually, this holds 
only for d > Ap + 1 = 5, while for d = 5a slightly different limit has to be taken 
switching off also the F{4) term. More precisely, we have to rescale r — > r/a and 
simultaneously b —>■ a 2 b. The resulting flat space solution is nothing but the 4 
dimensional instanton given by W = (1 — |6r 2 )/(l + |6r 2 ) considered as a static 
soliton of the d — 5 theory. Exactly the same function yields 'instanton' solutions 
for the F(2p) model in d — 4p + 1. 

The other endpoints of the curves b(d, a), obtained for large values of b, are 
more difficult to explain. These are different for d = 5 [S] and d = 6,7,8 
As was already observed in [5] for d = 5 there seems to be a critical value 
a c ~ 0.6 and there are obviously other turning points of b(d, a) besides the 
maximum a max ~ 0.75. For large values of b the solutions show a very distinctive 
behaviour. N decreases very quickly to some value iVo w 0.82, W stays close 
to one, while T 2 /r and U 2 /r tend to finite limits. We claim that the solution 
comes closer and closer to the conical f.p. (|HU|) as b — > oo. In Section El we gave 
an analytical description of these solutions leading to an asymptotic formula for 
the parameters a, c, and m for large b and to an independent determination of 
a c . The value a c = 0.595965 obtained this way agrees up to this precision with 
the value obtained with a regular solution for b = 1000. Employing Eq. (J30|) for 



16 



d = 5 we get iVo = y 2/3 ~ 0.8165 well in accordance with our numerical result 
(see Fig. |1J). Also the limiting value of y 1 / 4 = (1 — W 2 )/r 3 / 2 read of from Fig. |2] 
agrees well with the analytical result y)j^ = 2/ ^/3a c ~ 1.496. Fig. 01 displays the 
behaviour of b 3 (a c — a(b)) etc. for large values of b, which is in accordance with 
our derivation at the end of Section El As also mentioned in Section El the f.p. 
ceases to be oscillatory at d ~ 5.55, which is in accordance with the behaviour 
shown in Fig. 0] for d = 5.56. 

The f.p. exists only for d < d 2 with d 2 = (71 — V673)/8 ~ 5.632 given by 
Eq. (|32j) . Thus this f.p can only be relevant for the limiting solution up to this 
value, again well in accordance with our numerical results as shown in Fig. El 

For d > d 2 a different limiting behaviour manifests itself. As far as the 
numerical results for d = 6 are concerned, the function N seems to develop a 
double zero at some finite value of r, where the variable r increases without 
limit. This type of critical behaviour is well known from the limiting solution of 
the BM family, when the number of zeros goes to infinity and from gravitating 
monopoles as described in [18j . The critical solution runs into the RN f.p. and 
thus consists of a non-trivial (geodesically complete) interior solution and an 
exterior (geodesically incomplete) one with W — representing the exterior of 
the extremal RN black hole. 

Instead of looking for globally regular solutions one may as well look for 
solutions starting at r = and ending at the RN f.p. Fig. El shows corresponding 
curves b(d, a) as dashed lines. As discussed in Section El we have to rescale 
b — > bo?! 3 in order to obtain the BM solution of the F(A) model. The same 
rescaling has to be chosen for the solution running into the RN f.p. For d — 6 
there seems to be a bifurcation of the globally regular solutions with these RN 
ones at a k 0.015, although numerically it is clearly impossible to distinguish 
a genuine bifurcation from a finite but extremely small proximity. At any rate, 
it is rather obvious from Fig. El that the curve b(d, a) for d — 8 extends all the 
way to a = without any bifurcation with the RN curve, agreeing with the 
numerical results of 8J. Although for d = 7 this is not obvious from Fig. El 
our numerical results clearly exclude a bifurcation. In fact, numerically we can 
exclude a bifurcation down to d ~ 6.8. Thus one might expect some particular 
value d = df, for which the bifurcation occurs at a = and then for d < df, moves 
along some curve a(d) down to d = d 2 . In the following we shall give a convincing 
argument, that actually no such bifurcation happens for any d 2 < d < 9. This 
argument is based on an observation made for the BM solutions with many zeros, 
i.e. close to the limiting solution running into the RN f.p. There, solutions 
with n ^> 1 zeros show three characteristic regions. Region I extends from r = 
to r = 1, where the solutions develop more and more zeros of W and come close 
to the oscillating limiting solution running into the RN f.p., for which W decays 
like e~ T l 2 for r — > oo. In region II, between r = 1 and r = r n 1 the functions 
W and U first continue to decay, but then they start to grow like e T//2 , however, 
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still staying very small. In this interval the metric is well approximated by the 
extremal RN solution growing from N^O&tr = ltoNml for r>l. Thus 
the solution is well approximated solving the linearised equation for W and U 
in the background of the extremal RN solution. Finally, in region III between 
r = r n and r = oo the solutions stay close to the oscillating flat space solution 
for W running from r = oo into the f.p. with U = W = and N — k — 1. The 
relevant difference between the RN f.p. of the (U, W) system with N = and the 
flat one with N — 1 manifests itself in the real part of the eigenvalue being —1/2 
for iV = and +1/2 for N = 1. 

In the present case, if we assume a bifurcation with the solution running 
into the RN f.p. the same structure with three regions must form. Now we will 
argue that this leads to a contradiction. For a ^ the flat f.p. with W = is 
dominated by the p = 1 term at large r and thus is oscillatory in the interval 
\d-7\ < 2y/3 « 3.46. Since we are considering a family of solutions with one 
zero they cannot come close to this f.p. in the considered region with d « 6. 
This means that the curve b(d, a) has to turn back all the way to a = 0. There 
still remains the possibility that the bifurcation happens at a = for some or 
all d < db- Since b — > oo as a — > it is the F(4) theory that is relevant in 
this case. As found in Section H] the eigenvalues at the flat f.p. with W = are 
real only if d < 11 — 2^/5 ~ 6.528, thus d must obey this bound. Suppose there 
is some maximal value d m < 6.528 for which the p = 2 solution runs into the 
RN f.p., then it must come arbitrarily close to it for slightly larger values of d. 
As described above, the function W should then stay small in a large r-interval, 
while iV runs from iV w to iV w 1. From Eqs. (j25|26)l it follows that in the 
considered domain for d there are always two convergent modes in the (W, U) 
subsystem corresponding to two negative eigenvalues for the RN f.p. (r — > +oo) 
and two positive ones for the flat f.p. (r — ► — oo). However for solutions coming 
very close to both f.p.s the modes with the eigenvalues of larger modulus are 
suppressed, i.e. near the RN f.p. it is A+ N that dominates, while it is A^. at at the 
flat f.p. The best way to visualise the corresponding change of the (W, U) system 
moving from one f.p. to the other is to consider the quotient r] = rU/W obeying 
a Riccati equation in the linear approximation 

^ = 5_d + ( (1Q -d)N-K)ri-ri 2 , (48) 



with iV and k given by the extremal RN solution 

» T A, 2M rf-5 n| . , w 3 ( d-h \ 
iV =( 1 -^ + 4(9^) 3 - mth M =(93d)brf33)) ■ m 



d-3 
6 



and 

_ d - 3 + (5 - d)N 2 - 



(50) 



2N 

Supposing there is a solution with the correct behaviour at both f.p.s we should 
find a solution of the Riccati equation interpolating between A2 N at the RN f.p. 
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(r = To) and A^ at at the flat one (r — > oo). For d = 6 the Eq. (|48|) can be 
integrated in closed form resulting in a solution running from A+ N to A+ at instead 
to A^ at , excluding a bifurcation under the assumptions made. For other values 
of d in the relevant interval 5.632 < d < 6.528 numerical integration of Eq. (|4*%|) 
yields the same negative result. Therefore we conclude that the dashed and solid 
curves in Fig. El run back to different points at the b axis, although numerically 
we are unable to find regular solutions below d ~ 6.8. 

There remains the interval 4 < d < 5 for which no regular flat solutions exist 
according to the virial theorem. On the other hand there are self-gravitating 
solutions of the F(2) model generalising the BM solutions of d = 4. In fact there 
are again families b(d, a) of regular solutions tending to these BM type solutions 
in the limit 6 — > and a — > oo after suitable rescalings as discussed in Sectional 
Fig. |U1 shows some of the curves b(d, a). For b — > oo they tend to the conical f.p. 
as for 5 < d < d q . 

In order to further support our results for the 1-2-model we have also per- 
formed some numerical calculations for the 2-3- and the 2-4-models. Figs. [7| 
and |H1 show again agreement with our analytical results from Section El For 
d < Ap + 1 = 9 the curves b(d, a) have a BM type limit for a — > oo after suitable 
rescaling, while for d > 9 they start from a = with the flat space solutions with 
finite values of b. For d < ds ~ 8.63 resp. d < d& ~ 11.81 the curves approach a 
finite limit for a as b — > oo and the solutions run into the conical f.p., while for 
larger values of d the curves b(d, a) run back to a = approaching the BM type 
solutions of the F(6) resp. F{8) model after rescaling. 

The results on the f.p.s with W = and N = resp. N = 1 have another 
interesting consequence. Solutions of the p-g-model, which come close to the RN 
f.p. and for which iV — > 1 for large r come also close to the flat f.p. of the F(2p) 
model with W = (for r ^> 1 the F(2p) terms dominate those of the F(2q) part). 
Although the RN f.p. is non-oscillatory for p > 2, the flat one is oscillatory in 
the interval 4p + 3 — 2^/2p + 1 < d < Ap + 3 + 2^/2p + 1 with positive real part 
as long as d < 4p + l. Thus we expect to find regular solutions with any number 
of zeros of W accumulating at r — > oo for 4p + 3 — 2^/2p + 1 < d < Ap + 1. For 
p = 1 this interval is 3.54 < d < 5 (remember, however, that we assume d > 4), 
while for p = 2 we get 6.528 < d < 9. This is compatible with the observation of 
[S] that no multi-node solutions exist for the 1-2-model for d > 5. On the other 
hand for the 2-3-model multi-node solutions are expected to exist for the integer 
dimensions d = 7 and d — 8 and were actually found in our numerical analysis. 
Figs. 19*1 and ITU1 show solutions of the F(4) theory for d — 8 with two and three 
zeros of W. Solutions for d = 7 get much closer to the RN f.p. and are therefore 
more difficult to obtain numerically. Fig. ^2 shows a solution of the 2-3-model 
for d = 7 and a = 20 with two zeros of W. 

Since the eigenvalues of the f.p. have negative real part for d > 4p + 1 there 
are solutions with W — > for r — > oo (see Fig. fT2"jl . 
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7 Summary 



Solutions to a family of EYM models in higher dimensions have been studied. 
These models involve higher order YM curvature terms F(2p) characterised by 
integers p, p = 1 giving the standard YM theory. For reasons of simplicity we have 
restricted the number of such terms to two (p-q- models) depending on just one 
dimensionless parameter a. Studying the corresponding one-parameter families 
of (static, spherically symmetric) solutions allows us to capture the qualitative 
properties of EYM solutions in higher dimensions. We have restricted our study, 
for simplicity, to regular solutions only, knowing that corresponding black hole 
solutions exist. 

The simplest example of the family of theories we have considered (the 1-2- 
model) was recently studied in (HI E] devised specifically to suit spacetime dimen- 
sions 4 < d < 9. The more general p-g-models considered here suit also higher 
spacetime dimensions. The central part of the work here was aimed at under- 
standing qualitatively the results found from the numerical analysis of [HI El- This 
was achieved by analytic analysis supplemented with numerical studies including 
in addition to the previously considered case p = 1, q = 2 also the cases p — 2, 
q = 3 and p = 2, q = 4. The analytic part of the work was that of a fixed point 
analysis of nonlinear ordinary differential equations, employed previously [74^ IT5] 
for the familiar d — 4 BM solutions and magnetic monopoles. 

Like in the case of the gravitating monopoles |X2] of the EYMH system, one 
finds families of solutions parametrised by the dimensionless parameter a, which 
can be understood as a quotient of mass scales of the YM and the gravity parts 
of the theory. Since the a dependence of the solutions shows several distinctive 
behaviours for different dimensions it turned out to be particularly useful to 
vary the dimension parameter d continously, which is possible due to the explicit 
dependence of the radial field equations on d. This trick allowed us to relate the 
changes with d to changes in the fixed point structure of the equations. We find 
several different possibilities for the a dependence of the solutions that can be 
characterised by their behaviour for small b and for large b (see Figs. EJ El and|HJ). 
There are three possibilities for small b: 

• For 2q < d < Ap + 1 the parameter a decreases from infinity as b increases 
from zero. In the limit a — > oo the solutions approach a suitably rescaled 
BM type solution of the F(2p) model. 

• For d — Ap + 1, as b starts from zero the parameter a either increases from 
zero if Ap < 2q + 1, or decreases from infinity if Ap > 2q + 1. In the limit 
b —>■ the solutions approach the instanton of the flat F(2p) model. 

• For Ap + l<d<Aq+l the parameter a increases from zero as b starts 
from a finite value. For a — > the solutions approach a solution of the flat 
p-g-model. 
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And there are two possibilities for large b: 



• For 2q < d < d q the parameter a approaches a finite limiting value a c ^ 
as b — > oo. In the limit the solutions tend to a new type of fixed point, which 
we named 'conical' in view of the corresponding space-time singularity. 

• For d q < d < 4q + 1 the parameter a decreases to zero as b — > oo. For 
a — > the solutions approach the BM type solution of the F(2g) model. 
Note, however, that this limit cannot be explored numerically for values of 
d close to d q because the solutions come too close to the RN f.p. 

Clearly not all six combinations of these possibilities can be realised simultane- 
ously for one p-g-model, because some of the constraints on them are mutually 
exclusive. 

The existence of multi-node solutions is related to the existence of a f.p. with 
complex eigenvalues corresponding to a focal point. In the present case there 
is such a f.p. with vanishing YM potential W. Our analytical analysis reveals 
that for the p-g-models multi-node solutions are expected for dimensions in the 
interval 4p + 3 — 2y / 2p~+T < d < 4p + 1 perfectly compatible with our numerical 
results. 
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Figure 3: Rescaled parameters of the 1-2-model for d = 5 and large b 



25 



p=l q=2 d = 5.56 



0.4 - 



0.2 




(»- a nm) b3 - 73 - - (m-m lim )b 3 - 73 



_ Q g I I II | | | | | II | | | | | II | | | | | I L 

1 10 100 b iooo 

Figure 4: Rescaled parameters of the 1-2-model for d = 5.56 and large b 
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Figure 5: b as a function of a for regular solutions of the 1-2-model with d — 5, 
6, 7, and 8 
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Figure 8: b vs. a for regular solutions of the 2-4- model with different values of d 
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Figure 10: p — 2, d — 8 solution with three zeros of W 
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Figure 11: Solution of the 2-3-model for d = 7 and a = 20 with two zeros of W 
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Figure 12: Oscillating solutions of the p = 2 model for d = 8.9, 9, and 9.1 with 
decreasing, constant, and growing amplitude respectively 
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